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1.  EIGENSTRUCTURE  DECOMPOSITION  OF  TIME  SCALES 
IN  LINEAR-TIME  INVARIANT  SYSTEMS 


A.  Introduction 

Many  control  theory  concepts  are  valid  for  any  system  order,  however, 
their  actual  use  is  limited  to  low  order  models.  Large  scale  systems  result 
not  only  in  a  formidable  amount  of  computation,  but  also  in  ill-conditioned 
initial  and  two  point  boundary  value  problems.  The  interaction  of  fast  and 
slow  phenomena  in  high-order  systems  results  in  stiff  numerical  problems  which 
require  expensive  integration  routines.  The  singular  perturbations  approach 
to  decomposing  fast  and  slow  phenomena  involves  using  a  time-scale  separation 
technique.  In  this  case  a  reduced  order  "steady  state"  and  "boundary  layer" 
solutions  are  obtained  from  a  high  order  model.  Control  designs  and  simulations 
for  the  high  order  model  are  then  carried  out  on  the  reduced  order  subsystems. 

It  is  the  purpose  of  this  chapter  to  unify  and  extend  the  results 
of  previous  authors  [1-6]  and  attempt  to  provide  a  sense  of  completeness  to 
the  theory  of  time-scale  separation  in  linear  systems.  Our  concern  here  is 
the  linear  time  invariant  homogeneous  system 
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and  transform  it  into  either  the  form 


(1.1) 


(1.2) 


and/or  the  form 


2 


where 

inf  jo  (D^)  |  >  supja(Aoo)j 

(1.4) 

and 

inf |a(D*) |  >  sup | a (A*) | 

(1.5) 

and  any  such  system  (1.1)  which  has  this  property  is  said  to  satisfy  the  two- 
time-scale  property  for  dimensions  N  and  M. 

In  Section  B,  earlier  methods  of  time-scale  decomposition  are 
presented.  A  power  iteration  method  for  computing  the  dominant  left  eigen- 
space  of  a  matrix  is  used  to  show  the  equivalence  and  convergence  criteria 
of  the  past  schemes.  This  method  is  used  for  transforming  (1.1)  into  (1.2) 
satisfying  (1.4). 

Section  C  uses  a  power  iteration  for  computing  the  dominant  right 
eigenspace  of  a  matrix  to  show  the  existence  and  convergence  of  the  trans¬ 
formation  necessary  to  transform  (1.1)  into  (1.3)  satisfying  (1.5). 

Section  D  shows  the  duality  of  the  transformations  by  exploring 
the  eigenstructure  of  a  matrix  and  its  transpose.  , 

Section  E  completes  the  block  diagonalizations  of  (1.2)  and  (1.3) 
and  identifies  "fast"  and  "slow"  components  of  our  original  state  vectors. 

The  explicit  invertibility  of  our  transformation  matrices  is  shown.  This 
becomes  very  important  in  later  chapters. 

In  Section  F  we  consider  the  problem  of  properly  ordering  the  states. 
Finally,  in  Section  G,  we  give  an  illustrative  example. 


3 


B.  Earlier  Methods 

In  this  section  we  present  two  earlier  methods  for  transforming 
(1.1)  into  (1.2),  and  give  there  iterative  schemes, 
i)  Quasi  steady  state  method  [6] 

This  method  was  motivated  by  singularly  perturbed  models.  Here 
and  throughout  the  remainder,  it  is  assumed  that  the  states  are  ordered  and  N 
and  M  are  picked  such  that  D  ^  exists.  This  assumption  is  standard  in 
studies  of  singularly  perturbed  systems,  and,  as  will  become  apparent  later, 
this  assumption  is  not  restrictive. 

If  the  eigenvalues  of  D  are  such  that  the  real  parts  are  large  and 
negative,  then  the  homogeneous  solution  of  z  converges  to  a  steady-state 
rapidly.  If  this  convergence  is  assumed  to  be  instantaneous,  then  z  =  0  and 
this  quasi  steady  state  assumption  yields 

z  -  -D_1Cy  .  (1.6) 

s  s 

Next  we  try  to  remove  this  slow  part  of  z  by  introducing 

*  z  +  D  ^Cx  (1. 7) 

which  transforms  (1.1)  into 


(1.8) 


Repeating  steps  (1.6)  and  (1.7)  k  times  results  in  the  following 

\  ‘  Vl  *  \-lCk-ly  "o’2 


(1.9) 
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(1.10) 


where  the  subsystem  matrices  are  defined  as 


\  "  *k-l  "  BDk-lCk-l 

A  -  A 

0 

(1.11) 

ck  •  ClVA 

C  -  C 

0 

(1.12) 

-1 

(1.13) 

\  =  Vl  +  Dk-lCk-lB 

D-D. 

0 

Experimental  results,  motivated  by  singular  perturbations,  have 
converged  to  the  form  (1.2)  satisfying  spectral  property  (1.4). 
ii)  Algebraic  Riccati  equation  method  [3] 

In  [2, 3, 4, 7],  the  transformation  of  the  form 

n  *  z  +  Py  (1.14) 


is  proposed  in  an  attempt  to  transform  (1.1)  into  (1.2)  .  By  applying  (1.14) 
to  (1.1),  it  becomes 


A-BP 


B 


y 


C-DP  +  PA-PBP 


D  +  PB 


(1.15) 


The  problem  is  to  find  the  solution  P  to  the  Riccati  type  equation 

R(P)  =  C-DP  +  PA-PBP  =  0  (1.16) 

such  that  A-BP  and  D+PB  have  the  spectral  properties  (1.4).  Such  spectrum 
dependent  solutions  have  been  referred  to  as  "dichotomic"  [8].  We  will 
throughout  the  rest  of  this  report  continue  to  refer  to  this  solution  by  this 
label. 

Earlier  work  by  [4]  and  more  recent  work  by  [3]  have  resulted  in 


the  following  iterative  recursion  formula  for  obtaining  the  dichotomic 
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solution  to  (1.16) 


pk+i  ■  pk  + <D+pkB)  ■R(pk> 

P  -  D-1C 
o 


which  gives  the  subsystem  matrices  at  each  k  as 


(1.17) 


A  -  BP, 


Dk  "  D  +  Pk-lB 


A  “A 
o 


D-D 

o 


C.  -  R(P,  . )  -  C-DP,  ,  +  P,  .  A  -  P,  .BP,  ,  C  =  C. 
k  k-1  k-1  k-1  k-1  k-1  o 


(1.18) 


(1.19) 


(1.20) 


We  now  give  a  lemma  which  establishes  a  convergence  criterion  for 
(1.17).  In  the  process,  we  shall  show  that  (1.11)-(1.13)  and  (1.18)-(1.20) 
are  equivalent  at  every  iterate. 

Lemma  1:  Given  (1.1),  if  the  spectrum  is  concentrated  in  two  groups  of  M 
and  N  such  that 


supj  A^|  inf  I  A . | . 
i=l,N  j=l ,M 


(1.21) 


Then  under  mild  restrictions  on  the  initial  iterate  Pq,  (1.17)  will  converge 
to  the  dichotomic  solution  of  (1.16)  at  a  convergence  rate  of  e  ,  where 


suplxj  iei.N 
inf  |  Xj"|'  j€  1,M  ’ 


(1.22) 


Proof:  The  well  known  power  iteration  method  [9,10]  for  computing  a  M 
dimensional  basis  for  the  dominant  left  eigenspace  of  (1.1)  is  of  the  form 


fMk  -V’W  Nk-i5 


A  Bl 


(1.23) 


As  in  Ref. 
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R^  is  a  nonsingular  mxm  scaling  matrix  used,  for  example,  to  keep 

the  rows  of  [M^  N^]  strongly  independent  and  the  individual  components 

within  a  practical  range  of  computation  [9].  Many  methods  have  been  proposed 

for  selecting  the  sequence  of  R^'s  and  the  interested  reader  is  referred  to 

[9]  and  [11].  The  analytical  convergence  of  (1.23),  however,  is  independent 

of  R_  .  Thus,  under  condition  (1.21)  and  mild  conditions  on  [M  N  ] ,  it  is 
k  o  o 

known  that  (1.23)  converges  to  the  dominant  m-dimensional  left  eigenspace  of 

(1.1). 

Expressing  the  common  iterates  as 


\  ■  \-lA  +  Nk-lC  Mo  •  C  (1-24> 

"i;1  -  <d+\-iVib)'X-i  C*0'1-  a-25> 

Expressing  (1.24)  and  (1.25)  as  a  product, 

"kSc-  ®+<>k-il,>'1(c+VllMk-i«  (1-26) 

letting  Lk=N~1Mk  gives  (1.26)  as 

1^  *  (D  + Lk_1B)_1(C  +  Lk_1A)  Lq  -  D_1C  (1.27) 

which  is  equivalent  to  (1.17)  Vk>0. 

To  show  the  solution  is  dichotomic,  P  is  of  the  form 

P  =  N_1M.  (1.28) 

Without  loss  of  generality  we  can  take  [M  N]  =  [V^  V£]  whre  [V^  V^] 
are  the  M  left  generalized  eigenvectors  corresponding  to  the  M  dominant 
eigenvalues ,  thus 


TT’!' 


and. 


or, 


[Vl  V 


'a 

B 

r  j.  oi 

c 

D 

°  H 

> 

L  J 

41 

However,  from  (1.16) 


V  +  v2c  -  J2v1 


V  B  +  V2D  +  J2V2. 


tvi  v2] 


C  +  V21v1A  =  DV21v1  +  V21v1BV21vi 


which  leaves  both  (1.30)  and  (1.31)  in  the  form 


(1.29) 

(1.30) 

(1.31) 


V2(D  +  V2lviB)V21  =  J2 

verifying  the  desired  spectral  decomposition 

Corollary  2:  The  matrix  iterations  (1.11)-(1.13)  and  (1.18)-(1.20)  are 
equivalent  at  every  k.  Thus,  the  just  proved  convergence  properties  of 
(1.18)-(1.20)  are  propagated  to  (1.11)-(1.13) . 

Proof:  Substitution  of  (1.17)  into  (1.18)-(1.20)  results  in 

A-BPk_i  -  A-BPk.2-B(D  +  Pk_2B)-1-R(Pk.2) 

D  +  Pk-1B  ‘  D  +  Pk-2B  +  (0  +  Pk.2B)'1'R(Pk_2)-B 

.-1 


(1.32) 


Letting 


Ck  -  R(Pk_1)  =  (D  +  Pk_2B)  •R(Pk_2)-(A-BPk_1). 


\  =  A‘BPk-l 


D  +  Pk-1B 


\  =  Ck 


(1.33) 

(1.34) 

(1.35) 

(1.36) 

(1.37) 

(1.38) 


(1.33)  -  (1.35)  become 


a  * 

a,  ,  ~  By.  \  a.  , 

a  =  A 

(1.39 

k 

k-1  k-1  k-1 

-1 

0 

Yk  " 

Vi  +  Yk-iViB 

y  *  D 

0 

(1.40) 

ak  = 

Yk-lak-l\ 

a  =  C 

0 

(1.41) 

which  are  equivalent  to  (1.11)- (1.13)  Vk>0.  This  completes  the  proof. 

C.  Dual  Transformations 

Transforming  (1.1)  into  (1.2)  can  be  physically  interpreted  as 
removing  the  "slow"  components  of  the  z  states.  This  physical  interpretation 
motivated  the  Quassi-steady  state  iterations  of  the  previous  section. 

The  dual  to  this  procedure  involves  removing  the  fast  parts  of  the 
y  states.  Such  a  procedure  would  transform  (1.1)  into  (1.3)  satisfying 
condition  (1.5).  [6]  proposed  this  dual  procedure  which  lead  to  matrix 

recursions 


\  ■  Vi  -  \-iDk-ic 

A  =  A 

0 

(1.42) 

-1 

(1.43) 

Bk  =  \Bk-lDk-l 

B  =  B 

0 

Dk  =  Dk-1  +  CBk-lDk-l 

D  -  D. 

0 

(1.44) 

Likewise,  [4,5]  have  proposed  the  dual  to  the  Riccati  mathod  via  the  trans¬ 
formation 

C  ■  y  -  Pz  (1.45) 

which  leads  to  the  matrix  equations 


Ak"  A-Pk-1G 


Dk  “  D  +  CPk-l 


Bk-s(IW  - B  -  Pk-1D  +  ^k-i  -  pk-icpk-i 


A  =  A  (1.46) 
o 


D  =  D  (1.47) 
o 


B  =  B 
o 


(1.48) 
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where  ^  ~  ^  l 

Vi-pk  +  s<V(D+cV'  (1-49) 

P  =>  BD_1. 
o 

We  now  cite  a  lemma  [3]  which  is  the  dual  to  Lemma  1  and  establishes 

the  conditions  of  convergence  of  (1.49)  to  the  dichotomic  solution. 

Lemma  3:  If  the  spectrum  of  (1.1)  satisfies  (1.21),  then  under  mild 

restrictions  on  D  [3],  (1.49)  will  converge  to  the  dichotomic  solution  at  a 
o 

convergence  rate  of  (1.22). 

Proof:  The  well  known  power  iteration  method  [9,10]  for  computing  an  M 
dimensional  basis  for  the  dominant  right  eigenspace  of  (1.1)  is  of  the  form 


(1.50) 


serves  the  same  purpose  as  explained  in  Lemma  1. 
'Ml 

Thus,  under  mild  restrictions  on  N°  ,  it  is  well  known  that  (1.50)  converges 

o_ 

to  the  dominant  M-dimensional  right  eigenspace  of  (1.1). 

Expressing  (1.50)  as 


where  the  scaling  matrix 


Vl  •  <B  +  AV«>+CPk) 

p  -  bd"1 
0 


(1.51) 


(1.52) 


(1.53) 
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which  is  equivalent  to  (1.49)  ^>0.  Proving  P  is  dichotomic  is  carried  out 

as  in  Lemma  1  and  can  be  seen  in  [3,4]. 

Corollary  4:  The  matrix  recursions  (1.42)-(1.44)  and  (1.46)— (1.48)  are 


equivalent 

at  every  k,  thus,  the  convergence  properties  of  (1.52) 

are  propa- 

gated  to  (1 

•42)-(1.44) . 

Proof:  Substitution  of  (1.49)  into  (1.46)-(1.48)  gives 

A-Pk-lC  -  4-Pk-2C-S<Pk-2)(D  +  CPk-2)'lc 

(1.54) 

0+aW  •  °  +  CPk-2  +  CS(\-2>(I>+CPlt-2>'1 

(1.55) 

Bfc  -  S(P][.l)  =  (A-Pk.1C)S(Pk.2)(D-CPt.2)-1. 

(1.56) 

Letting 

\  =  A  ~  Pk-lC 

(1.57) 

Yk  *  D  +  CPk,L 

(1.58) 

ek*Bk 

(1.59) 

gives 

\  =  vr8k-iVic  ao " A 

(1.60) 

Yk  =  Vi  +  cek-iVi  Yo  * D 

(1.61) 

’  Vk-iVi  so  ■ B 

(1.62) 

which  are  equivalent  to  (1.42)-(1.44)  . 


D.  Some  Comments  on  the  Eigenstructure  of  a  Matrix  and  its  Transpose 
Consider  the  two  system 

x  *  Fx  (1.63) 


T»  a 

x  -  Fx. 


(1.64) 
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of 

(1 

.1), 

(1. 

63) 

*y 

"a 

B  ' 

y  ’ 

z  . 

.c 

D  . 

_  z 

J  - 

r.T 

T  1 

r  i 

y 

A 

c 

y 

s 

_T 

_z  . 

LB 

D  J 

z 

(1.65) 


(1.66) 


Now,  let  us  try  to  decompose  time  scales  in  (1.65)  and  (1.66) 
independently.  First,  in  (1.65),  let 

w  »  z  +  Ly 

which  transforms  (1.65)  into 


w 


(A-BL)  B 

C-DL  +  LA-LBL  D  +  LB| 


where  we  seek  L  as  the  dichotomic  solution  to 

R(L)  *  C-DL  +  LA  -  LBL  =  0 . 

Now,  consider  the  transformation 

x  *  y  -  Lz. 

applied  to  (1.66)  which  gives 


at-lbt  atl  +  ct-lbtl-ldt 

T  T  T 

B  D  +B  L 


(1.67) 


(1.68) 


(1.69) 


where  we  seek  L  as  the  dichotomic  solution  to 


R(L)  =■  CT-LDT  +  ATL-LBTL  =  0. 


(1.70) 


One  method  of  obtaining  a  solution  to  (1.70)  is  through  the  power  iteration 


Vii 


rAT 


,T-, 


T  T 
B  D 


L\ 


(1.71) 
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Letting  the  dichotomic  solution  of  (1.70)  is  the  equilibrium 

solution  of 

Lk+1  “  (A^+C^CB^+D1)"1  (1.72) 


T  T' 
L  =  CD1 
o 


-1 


To  obtain  the  dichotomic  solution  for  (1.68)  first  note  that 


RT<t>  .(?-tVtAT-tW.o. 


(1.73) 


Thus,  the  solution  to  (1.68)  obtained  using  (1.73)  is  the  transpose  of  the 
dichotomic  solution  to  (1.70). 


\  ‘  (MkA1)T 

T  i  T 

■K  K 

A 

L  -DC 


and  the  power  iteration  (1.71)  becomes 


or, 


^\+l  Nk+1^ 


[Mk 


(1.74) 


(1.75) 


(1.76) 


(1.77) 


The  conclusions  are: 

i)  The  power  iterations  for  computing  the  right  dominant  eigenspace  for 
(1.66)  are  equivalent  to  computing  the  left  dominant  digenspace  of  (1.65) 

ii)  Removing  the  "slow"  from  z  is  equivalent  to  removing  the  "fast"  from  y. 

iii)  L-LT. 


iv)  The  subsystem  matrices  obtained  from  the  spectrum  decomposition  are 
related  by 
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T  T  ~  T 

(A  -  LB  )  =  (A  -  BL) 

T  T  -  T 

(D  +BL)  =  (D  +  LB) 

Analogus  results  hold  when  the  transformations 

w  =*  z  +  Ly  ,  x  *  y  -  Lz 

are  applied  to  (1.65)  and  (1.66)  respectively. 


E.  Block  Diagonalization  and  Identification  of  Fast  and  Slow  State  Vector 
Components 

Once  we  have  transformed  (1.1)  into  (1.2)  or  (1.3)  satisfying 
conditions  (1.4)  or  (1.5)  respectively,  block  diagonalization  is  always 
possible. 

Consider  form  (1.2) ,  and  the  transformation  (1.14)  used  to  obtain 
this  form.  The  dichotomic  solution  matrix  P  is  of  the  form 

P  =N_1M 


where  the  rows  of  [M  N]  span  the  dominant  left  eigenspace  of  (1.1).  Thus, 
the  exact  form  of  (2)  is 


_  — 

-i 

—  — 

y 

3> 

A  -  BP 

B 

y 

w 

0 

D  +  PB_ 

w 

(1.78) 


Now,  let  x*y-Qw.  This  leaves  (1.78)  in  the  form 


X 

A  -  BP 

(A  -  BP)Q  -Q(D  +  PB)  +B 

X 

• 

w 

0 

D  +  PB  _ 

w 

(1.79) 


Thus,  we  seek  Q  to  satisfy  the  Lyapunov  type  equation 


(A  -  BP)Q  -Q(D+PB)  +B  =0 


(1.80) 


14 


where 
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x(0-e(A-BP)txo  xo  =  (I-QP)yo-QZo 

w(t)  *  e^  +  PB^Cw  w  *  Py  +  z 

'  7  q  ©  o  o 

Such  decompositions  will  become  more  important  when  we  consider  singularly 
perturbed  systems  in  the  next  chapter.  There,  the  fast  and  slow  components 
take  on  the  names  of  "Boundary  Layer"  and  "Steady  State  Components". 

Now  consider  form  (1.3)  and  the  transformation  (1.45)  used  to  obtain 

A 

this  form.  The  Dichotomic  solution  matrix  P  is  of  the  form 


P  •  MN~ 1  (1.87) 

T 

Where  the  columns  of  [M  N]  span  the  dominant  right  eigenspace  of 
(1).  Thus,  the  exact  form  of  (1.3)  is 

|  A  -  PC  0  x 

C  D  +  CP  z 


“  — 

X 

9 

z 

—  _J 

(1.88) 


Now,  let  w  =  z+Qx.  This  transforms  (1.88)  into  the  form 


X 

A  -  PC 

0 

X 

w 

Q(A  -  PC)  -  (D  +  CP)Q  +  C 

D+CP_ 

w 

(1.89) 


Thus  we  seek  Q  to  satisfy  the  Lyapunov  type  equation 

Q(A  -  PC)  -  (D  +  CD)Q  +  C  »  0  (1.90) 

such  a  Q  will  always  exist  since 

ct  (A  -  DC )  n  a  (D  +  CP)  *  0  (1.91) 


Again,  (1.90)  may  be  solved  iteratively  or  algebraically.  Applying  the 
dominant  left  eigenspace  iterations  to  (1.88)  convergence  is  assured.  This 
iteration  takes  the  form 
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Qk+i*  (D  +  CP)"1  ’  (C+Qk*(A  -  PC))  Qq  =D_lC 
The  resulting  system  is  of  the  form 


—  - 

X 

“a  -  PC 

0 

"y“ 

• 

w 

0 

D  +  CP_ 

z 

and  the  composite  transformation  is 


—  *— 

— 

**  - 

• 

X 

3 

I 

-p 

y 

w 

_Q 

I  -  QP_ 

z 

with  also  possesses  the  explicit  inverse 


y 

I  -  PQ 

P 

X 

z 

_  < 

w 

(1.92) 


(1.93) 


(1.94) 


(1.95) 


Thus,  we  have  again  decomposed  the  y  and  z  state  vectors  into 
"Fast"  and  "Slow"  components.  Namely 


y 

z 


(I  -  PQ)x  +  Pw 


-Qx+w 


2slow 


yslow  +  ^fast 
+  z£ast 


(1.96) 

(1.97) 


where 

x(t)=e^A  TC>x  x  =  y  -  Pz 
o  o  J  o  o 

w(t)=e^P  +  ^P^w  w  *Qy  +(I-QP)z 
o  o  yoNX/o 

The  relationships  between  various  fast  and  slow  components  of  (1.1)  will  be 
made  apparant  in  the  next  chapter. 


F.  Ordering  of  States  Variables 

In  both  the  dominant  left  and  right  eigenspace  iterations,  it  was 
assumed  that  if  (L.21)  held,  then 
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(D+PkB) 


(1.98) 

(D+CPk)'L  (1.99) 

exists  for  every  k.  (1.17)  is  used  to  obtain  the  dichotomic  solution  of 
(1.16).  This  solution  was  shown  to  be  of  the  form 

(1.100) 


P-V21V21 


Where  [1^  ^22^  are  eigenvectors  corresponding  to  the  dominant 
eigenvectors  of  (1.1).  If  V22  does  not  exist,  then  neither  will  (1.100).  To 
analyse  this  problem  look  at 


11 


'21 


'12 


-22 


A11  A12 
A2l  A22 


A. 


0 

A„ 


11 


'21 


'12 


22 


or, 


VUA11+V12A21=A1V11 


V11A12  + V12A22  * A1V12 


V21A11  +V22A21  ~  V^l 


V21A12+V22A22  S‘A2V22 


(1.101) 

(1.102) 

(1.103) 

(1.104) 


[14]  has  shown  that  for  arbitrary  N,  M,  there  will  exist  an  ordering  of  the 
states  (eigenvectors)  such  that  V22  exists.  However,  [14]  does  not  guarantee 
that  the  resulting  spectrums  of  and  satisfy  (1.21). 

In  the  next  chapter  we  will  consider  singularly  perturbed  systems. 
In  this  case,  the  system  matrix  possesses  certain  explicit  properties  that 
make  it  well  suited  for  (1.98)  and  (1.99).  To  see  this,  assume 

llA22li  “  llA2l11  “  llA2H  “  0(1/^)  (1. 105) 


Then  (1.104)  can  be  written  as 


r 
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V22A22+0(^>  “A2V22  (1.106) 

Thus,  V^2  is  0(p.)  close  to  being  the  left  eigenvector  matrix  for 
A£2«  Using  the  relations  (1.105),  (1.98)  becomes 

(D  +  ^PkB)_1  (1.107) 

Using  bounds  on  ij.  established  by  [4]  on  P,  we  can  establish 
sufficient  bounds  on  n  guaranteeing  the  existence  of  (1.108)  V  ,  thus 
guaranteeing  the  existence  of  ^22'  ^us,  c^e  ordering  of  states  such  as  to 
fit  a  singularly  perturbed  format  is  very  important  to  the  success  of  (1.97) 
and  (1.98).  [3]  gives  an  algorithm  for  putting  an  arbitrary  system  into 

singularly  perturbed  format  and  should  be  used  before  any  multiple  time 
scale  analysis  is  attempted. 

Likewise,  to  analyse  (1.98), 

A11  A12 
__A2l  A22 

or. 


using  (1.105),  (1.112) 

and  again  the  desired  properties  become  obvious. 

It  is  interesting  to  note  that  if  and  exist,  then  using 


E11  E12  E11 

33 

E21  E22  _ E2l 

E12  A1 
'2dh 

0 

A2_ 

(1.108) 

aUEU+a12E2i  =eual 

(1.109) 

A 1 1E 1 2  +A12E22  *  E12A2 

(1.110) 

A21E11  +A22E21  =E21A1 

(1.111) 

A21E12+A22E22  =E22A2 

(1.112) 

and  (1.110)  yield 

“*"^“22^22  =  ^*22^2 

(1.113) 

1 
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VE  =  I, 


(1.114) 


we  obtain 


V 

E 


22V21 

E’1 

12^22 


'E2lEll 


(1.U5) 

(1.116) 


The  quantities  on  the  left  of  (1.15)  and  (1.116)  are  the  equilibruim  solutions 
to  the  dominant  left  and  right  eigenspace  iterations  respectively.  Thus 
the  fact  that  we  ordered  our  states  so  that  there  were  M  dominant  eigenvalues 
and  the  z  state,  were  fast  was  arbitrary  in  that  by  reordering  our  states, 
the  right  eigenspace  method  may  be  used  for  (1.2)  and  the  left  for  (1.3). 

The  existence  of  all  of  these  inverses  simulataneously  is  unlikely.  How¬ 
ever,  the  relations  (1.115)  and  (1.116)  do  help  to  further  clarify  the  dual 
nature  of  our  left  and  right  eigenspace  transformations. 


G.  Example  -  Decomposition  of  States  into  Fast  and  Slow  Components 

In  [13],  the  8th  order  model  of  an  isolated  mixed  power  system  is 

given  as 


-.2 

0 

0 

0 

0 

0 

0 

0 

4.75 

-5 

0 

0 

0 

0 

0 

0 

0 

.16667 

-.16667 

0 

0 

0 

0 

o 

0 

0 

2 

-2 

0 

0 

0 

0 

0 

1 

o 

00 

-.07467 

-.112 

-3.9944 

10 

-.92778 

-9.1 

0 

0 

0 

0 

.2 

-.5 

0 

0 

0 

0 

0 

0 

0 

0 

-1.39 

-.278 

0 

.01 

.0093 

.014 

-.06319 

0 

.11597 

-.112361 

(1.117) 

r 


using 

a  permutation  of 

x  *  Py 

P  = 

(el,e3’e65 

W  V 

W 

gives 

— 

_ 

-.2 

0 

0 

0 

0 

0 

0 

0 

0 

-.167 

0 

0 

0 

0 

0 

.167 

0 

0 

-.5 

0 

.2 

0 

0 

0 

y  a 

0 

.009 

0 

-.112 

-.063 

.014 

.116 

01 

0 

-.075  10.00 

-9.101 

-3.994  - 

.112 

-.927 

-.08 

0 

2.00 

0 

0 

0 

-2.00 

0 

0 

0 

0 

0 

-.277 

1.319 

0 

-1.386 

0 

4.75 

0 

0 

0 

0 

0 

0 

-5.00 

(1.118) 

The  eigenvalues  of  (1.118)  are 


Using  N  ■  M  *  4, 


eigenspace  iterations  ve 


A  -  P  B  - 
o 


-0.20000 

0.15834 

-0.00312 

0.00877 


-1.3884147 

+ 

0. 0000000 J 

-0.1291288 

+ 

0.2124795J 

-0.1291288 

jm 

0.2124795J 

-4.3489879 

+ 

0.0000000J 

(1.119) 

-2.0000000 

+ 

■  O.OOOOOOOJ 

-0.1666700 

+ 

0. 0000000 J 

-5.0000000 

+ 

O.OOOOOOOJ 

-0.2000000 

+ 

O.OOOOOOOJ 

we  obtain  an 

e 

of  .1792.  Using  the 

dominant 

left 

obtain 

0.00000 

0.00000  0.00000 

-0. 16667 

0.00000  0.00000 

(1.120) 

-0,00766 

- 

0.08981  -0.36571 

0.02153  0.09635  -0.22145 


i 

i 


which  has  eigenvalues 
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-0.15563  +  0.17579J 
-0.15563  -  0.17579J 
-0.16667  +  0.00000J 
-0.20000  +  O.OOOOOJ 


D  +  BP  = 
o 


-4.52014 

0.00000 

0.80733 

0.00000 


-0.08640 

-2.00000 

0.02712 

0.00000 


-0.71572 

0.00000 

-1.16426 

0.00000 


-0.05533 

-0.16667 

0.02543 

-5.00000 


which  has  eigenvalues 


-4.33808  +  O.OOOOOJ 
-1.34632  +  O.OOOOOJ 
-2.00000  +  O.OOOOOJ 
-5.00000  +  O.OOOOOJ 


using  the  dominant  right  eigenspace  iterations  we  get 


-0.20000 

0.00000 

0.00000 

0.00000 

A -CP  = 

0.15834 

-0.16667 

0.00000 

0.00000 

o 

-0.00312 

-0.00766 

-0.08981 

-0.36571 

_0. 00877 

0.02153 

0.09635 

-0.22145_ 

which  has  eigenvalues 

-0.15563  + 

0.17579J 

-0.15563  - 

0.17579J 

-0.16667  + 

O.OOOOOJ 

-0.20000  + 

O.OOOOOJ 

-4.31691 

-0.03023 

0.04757 

-0.05415 

D  +  P  C  ■ 

0.00000 

-2.00000 

0.00000 

-0.06667 

0 

1.32208 

0.00179 

-1.36749 

0.00051 

0.00000 

0.00000 

0.00000 

-5.00000 

(1.121) 


(1.122) 


(1.123) 


(1.124) 


d.125) 


D  +  P  C  = 
o 


(1.126) 
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Which  has  eigenvalues 


-4.33808  +  0.00000J 
-1.34632  +  0.00000J 
-2.00000  +  O.OOOOOJ 
-5.00000  +  O.OOOOOJ 

To  show  how  this  accuracy  may  be  improved,  after  two  iterations  we 


-0.20000 

0.00000 

0.00000 

o.ooooo" 

A  -  PB  = 

0.16492 

-0.16667 

0.00000 

0.00000 

-0.00233 

-0.00741 

-0.08420 

-0.36234 

0.00979 

0.02836 

0.13210 

-0.17633 

with 

eigenvalues 

-0.13027  + 

0.21388J 

-0.13027  - 

0.21388J 

-0.16667  + 

O.OOOOOJ 

-0.20000  + 

O.OOOOOJ 

-4.52468 

-0.08664 

-0.71767 

-0.0557~ 

D  +  BP  = 

0.00000 

-2.00000 

0.00000 

-0.18172 

0.76779 

0.02154 

-1.21045 

0.01334 

_0. 00000 

0.00000 

0.00000 

-5.00000 

with 

eigenvalues 

-4.34912  + 

O.OOOOOJ 

-1.38601  + 

O.OOOOOJ 

-2.00000  + 

O.OOOOOJ 

-5.00000  + 

O.OOOOOJ 

— 

_ 

-0. 2QCCQ 

0.00000 

0.00000 

0.00000 

A  -  CP  - 

0.16379 

-0.16667 

0.00000 

0.00000 

-0.00237 

-0.00609 

0.00444 

-0.45814 

0.00821 

0.02283 

0.13946 

-0.26497 

(1.127) 


obtain 


(1.128) 


(1.129) 


(1.130) 


(1.131) 


(1.132) 
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with  eigenvalues 


-0.13027 

+ 

0.21388J 

i 

-0.13027 

- 

0.21388J 

(1.133) 

-0.16667 

+ 

0. 00000 J 

-0.20000 

+ 

O.OOOOOJ 

j 

i 

-4.37192 

-0.03420 

-0.05145 

-0.05565 

j 

\ 

A 

D  +  PC  = 

0.00000 

-2 . 00000 

0.00000 

-0.06896 

(1.134) 

! 

i 

1.32327 

0.00202 

-1.36321 

0.00048 

\ 

0.00000 

0.00000 

0.00000 

-5.00000 

1 

with  eigenvalues 

\ 

\ 

1 

i 

-4.34912 

+ 

O.OOOOOJ 

j 

i 

\ 

-1.38601 

+ 

O.OOOOOJ 

(1.135) 

j 

i 

1 

j 

-2.00000 

+ 

O.OOOOOJ 

-5.00000 

+ 

O.OOOOOJ 

\ 

! 

We  now  give  graph  of  the  various  states  along  with  their  fast  and 
slow  components  using  (1.85),  (1.86),  (1.96)  and  (1.97)  for  both  the  left 
and  right  eigenspace  decompositions.  The  plots  will  be  based  on  the  PQ  and 

A 

Pq  iterates.  On  the  graphs  of  the  individual  components,  the  following 
legend  will  be  in  effect 

ACTUAL  STATE  - 

SLOW  COMPONENT  . 

FAST  COMPONENT . 

On  the  graphs  of  the  actual  state  versus  the  approximated  state 

ACTUAL  STATE  - 

APPROXIMATED  STATE  . 

The  plots  appear  on  the  next  several  pages.  The  system  is  perturbed 


with  an  initial  state  vector  of 


:seconds> 


Figure  1.2.  State  1  and  added  components  using  right  eigenspace  iterations. 


CSCN003S : 


TIME  C  SECONDS-}/ 

Figure  1.6.  State  2  and  added  components  using  right  eigenspace  iterations. 


T  CO  CVJ  *-  cs 

I 


WH  <HUJ  <0  •  •  •  Cx  _llii 

_  j 


:  SECONDS > 


State  3  and  added  components  using  left  eigenspace  iterations. 
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<0 


QKbJ 


Figure  1.9.  State  3  and  componen 


<0  H  <  H-  UJ  CO 


q  au 


Figure  1.10.  State  3  and  added  components  using  right  eigenspace  iterations 
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r\ 

V) 

a 

z 

o 

o 

Li 

V) 

‘w* 

LJ 

Z 

H 

I- 


<0  H  <  H  LxJ  '? 


Cx  -ltd 


Figure  1.11.  State  4  and  components  using  left  eigenspace  iterations. 


Figure  1.12.  State  4  and  added  components  using  left  elgenspace  iterations. 
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I  SECONDS} 


rations . 


1/5 


Figure  1.17.  State  5  and  component 


Figure  1.19.  State  6  and  compone 


Figure  1.20.  State  6  and  added  components  using  left  eigenspace  iterations. 


I 
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CS0N033S: 


Figure  1.23.  State  7  and  components  using  left  eigenspace  iterations. 
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<0  H  <  h-  lu  r> 


Q  _i  UJ 


Figure  1.24.  State  7  and  added 


Figure  1.29.  State  8  and 


Figure  1.30.  State  8  and  added 
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II.  ASYMPTOTIC  SERIES  DECOMPOSITION  OF  TIME-SCALES 
IN  LINEAR  TIME- INVARIANT  SYSTEMS 

A.  Introduction 

When  small  parameters  are  present  in  differential  equations  defining 
initial  or  boundary  value  problems,  one  of  the  popular  methods  of  solution 
[15]  is  to  obtain  an  asymptotic  power  series  expansion  of  the  solution.  Such 
techniques  have  been  well  documented  and  can  produce  approximate  solutions 
to  problems  where  otherwise  analytic  explicit  solution  are  impossible  or 
exact  numerical  solutions  are  computationally  not  practical.  Such  systems 
are  of  the  form 

*  -  f (X,t,e)x(0)  =  Xq  (2.1) 

and  a  solution  is  sought  of  the  form 

X(t)  =X°(t)  +  Sx1(t)  + .  (2.2) 

when  such  an  expansion  converges  uniformly  in  X  as  e-*  0  we  have  a  regular 
perturbation  problem  [16],  If  there  is  a  region  of  nonuniformity,  usually 
at  one  of  the  boundaries,  we  have  a  singular  perturbation  problem.  In  most 
cases,  the  dynamics  of  the  solution  vector  within  this  region  of  nonuniform 
convergence  involve  fast  transients  or  the  so  called  "Boundary  Layer 
Phenomena".  Thus,  such  singularly  perturbed  systems  [16]  are  said  to  possess 
an  inherent  two- time- scale  property  characterized  by  a  steady  state  or  outer 
solution  which  is  defined  by  the  regions  of  uniform  convergence  of  (2.2), 
and  the  boundary  layer  or  inner  solution  where  a  stretched  time  variable  is 
usually  introduced  in  order  to  achieve  convergence  on  the  total  time  interval. 
In  the  linear  case  such  systems  take  the  form 
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y  =  Ay  +  Bz  y(t)  =  yQ 
p>z  ■  Cy  +  Dz  z(Cq)  =  zq 


(2.3) 


Much  work  has  been  done  to  exploit  the  multiple  time  scale 
property  of  (1.1)  when  designing  regulators,  pole  placement,  reduced  order 
modeling,  etc.  [17,18,19].  As  a  result,  for  a  system  which  is  known  to 
have  fast  and  slow  phenomena,  the  systems  engineer  is  motivated  to  permute 
the  state  in  order  to  attain  the  above  structure  and  take  advantage  of  these 
decomposition  techniques.  It  is  the  purpose  of  this  chapter  to  use  multiple 
time  scale  asymptotic  expansions  to  obtain  a  "Steady  State"  and  "Boundary 
Layer"  decomposition  in  (2.3)  and  see  how  this  decomposition  compares  to  our 
eigenstructure  decompositions  of  Chapter  1. 

In  Section  B  we  obtain  power  series  representations  of  our  dichotomic 
transformation  matrices  P  and  P.  First  an  asymptotic  power  series  is  derived, 
then  its  equivalence  to  a  convergent  MacLauren  series  is  established. 

In  Section  C  we  derive  important  relationships  between  various 
fundamental  sets  of  solution  of  (2.3)  and  any  system  satisfying  (1.21).  These 
fundamental  sets  are  based  on  our  reduced  order  subsystem  matrices  and  the 

*  /S 

dichotomic  transformation  matrices  P,  P,  Q,  and  Q. 

In  Section  D  we  use  Vasil 'evas  method  of  matched  asymptotic  expansions 
to  obtain  the  "Boundary  Layer"  and  "Steady  State"  components  of  the  solution 
vectors  y(t)  and  z(t).  It  is  shown  that  this  decomposition  is  equivalent 
to  the  eigenstructure  decompositions  of  Chapter  1  by  using  one  of  the  funda¬ 
mental  matrices  established  in  Section  C. 

Section  E  discusses  some  computational  simplifications  to  the 
dominant  left  and  right  eigenspace  iterations  based  on  system  (2.3).  The 
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simplifications  involve  eliminating  the  necessity  to  take  an  inverse  at 
every  iteration.  The  simplified  iterations  were  originally  proposed  by  [4] 
where  a  contraction  mapping  argument  was  used  to  obtain  bounds  for  convert 
gence.  Using  a  different  approach,  less  conservative  bounds  are  obtained. 

Finally,  in  Section  F  we  give  an  example  which  highlights  many  of 
the  important  results  of  this  Chapter. 
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B.  Series  Solutions  to  Ricatti  Iterations 

For  proper  spectral  decomposition  and  dimensions  m  and  n  in  (1), 
it  was  shown  in  section  C  of  Chapter  1,  that  the  following  matrix 
recursion  equation 

fk+i  *  pk  -  <D+Pt«'1-  <DVc+Vpk-pkA)  <2-4> 

P  =  D_1C 
o 

will  converge  to  the  dichotomic  solution  of 

R(P)  =  DP  -  C +PBP  -  PA  =  0  (2.5) 

when  (1.1)  is  in  the  form  of  a  singularly  perturbed  model,  (2.4) 
becomes 

Pk+1  =  Pk  '  (D  +  aPkB) " 1  '  (DPk  ‘  C  +  uPkBPk  *  “PkA)  (2  • 6  > 

P  =  D_1C 
o 

In  a  later  section  in  this  chapter,  we  will  use  the  method  of 
matched  asymptotic  expansions  in  an  attempt  to  decompose  (2.3)  into 
"steady-state"  and  "boundary  layer"  subsystems.  It  will  be  necessary  in 
this  derivation  to  obtain  a  power  series  solution  to  (2.5)  in  the  form 

P  =  P°  +P.P1  +m.2P2  + .  (2.7) 

From  (2.6),  it  is  important  to  note  that  such  an  expansion 
converges  as  m>  -*  0  uniformly  in  K.  Thus,  we  have  a  regular  perturbation 
problem  [16]. 

Expressing  (2.6)  in  a  more  convienient  form 


<D+u,PkWPk+1 


C  +  UP,  A 
k 


(2.8) 
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We  now  substitute  in  the  formal  series  (2.7)  and  obtain 
[(D+u(P°+^pJ  +  ....)B)(P®+1  +  up£+1 +  ....)] 

-  C  +  u(P°+Mj£  +  ..)A 

Equating  like  powers  of  u  we  obtain  the  so  called  "equations 
of  the  variations"  [20] . 


(2.9) 

(2.10) 

(2.11) 


(2.12) 


(2.13) 


(2.14) 


Vk  S  0 


(2.15) 


Likewise,  the  solution  to  (2.10)  is  found  from 


4i  ■  •  D'lpX+i +  °'lpkA 


(2.16) 


However,  since  P°  =  P.°  Vk  s  0 
k+1  k 


■  D  ^P°BP°  +  D  ^P°A  =  constant  matrix 


Thus , 


P,  =  0  k  »  0 
k 


-loo  “1  o 

*  D  PkBPk  +  D  PfcA  Vk  1 

Continuing  in  this  manner,  the  equilibrium  solution  for  the  N1-*1 
variation  becomes 


(2.17) 


=  -  D_1  £  P^‘1_jBP^  +  D_1P^_1A 

k  .  „  k  k  k 

j=0 


(2.18) 


Thus,  we  have  a  well  defined  method  of  generating  all  terms  in  the  series. 
However,  to  prove  that  the  series  (2.7)  does  asymptotically  solve  the 
matrix  recursion  (2.8)  as  \i  -»  0,  let 


P(H)  =  2  p-V  +  R(u)  N  £  0 

j-o 


(2.19) 


N+l 

We  shall  show  that  the  remainder  term  is  unique  and  0(u  ). 

Substituting  (2.19)  into  (2.8) 

D(pk+i+upi+i  + . +“t,p£+i  +  iw>‘» 

+  “<pk+upk+--  +“i,<+Rk(u>)B<pk+1+“pi+i+--+“'‘pk+i+\«  <“> 


u(P°+uPk  +  M.2Pk  +  .  .  +y,NPk+Rk(M.))A  +C 


(2.20) 


6L 


Using  the  series  solutions  to  P^,  j  =  1,  N,  (2.20)  reduces  to 


DRk+i(^  +“  B  +pk(ti)BW 


+  +  V‘‘)BPiU> 


+Rk<">B<+l> 


+  uN+1f(p°,pJ,  ....  i£)  =  m^cujA 


(2.21) 


Collecting  terms , 


R^fo)  =  (D+u(P°+g,P^  + 


N  N  -1 

+  4  p£)B  +  liR^)  . 


'“\<u)B(pt« +“piU + •  •  • + “X+i5 + “V**' 


+  0(,,i+1) 


However,  since  Rq  =  0,  (2.22)  can  be  uniquely  solved  recursively 


i,e’  Rj_  =  (D  +y.P°B)_1  -  0  +  0(nN+1)  =  0(uN+1) 


R  =  (D+u  I  B  +  0(  uN+2))_1  0(4N+2) 
j=0 


^  N+l,  N+l, 

+  0(u  )  =  0(4  ) 


and  the  unique  solution  is  given  by 


(2.22) 


* 
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R(u)  =  In  R^yt)  =  0(nN+1) 

1C-® 


Thus,  the  solution  to  (2.22)  can  be  expressed  as 


IN  .  .  „ 

P(U)  =  T.  uV  +0(u  ) 
1=0 

VN  2  0 


or  using  the  notation  of  [15], 


(2.23) 


(2.24) 


P(U)  ~  £  u.jPj  (2.25) 

j=0 

Asymptotic  power  serie  like  (2.25),  in  general  do  not  converge 
for  any  fixed  nonzero  value  of  ij,.  If  however,  we  can  show  the  function 
P^Oi)  to  be  holomorphic  in  y,  in  some  neighborhood  of  the  origin  in  the 
complex  y  plane,  then  we  know  that  for  sufficiently  small  y,  P(y) 
possesses  a  convergent  MacLauren  series  expansion.  And  by  the  uniqueness 
of  this  series,  it  is  equivalent  to  (2.25)  [15]. 


Lemma  5 


Let 

If  F  satisfies 
the  orgin  in  y 


F 


A 


cNXN 


and 


y  e  c . 


|_C  DJ 

property  (1.21),  then  there  exists  a  neighborhood  S  of 
space  such  that  P^(y)  is  holomorphic  in  S  Vk  s  0. 
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Proof : 


To  prove  analyticity  of  P  in  a  define  region  S  of  the  u  plane, 

rC 


it  is  both  necessary  and  sufficient  to  show 


exists,  and  is  continuous  in  S 


We  show  this  by  induction 


P  *  D  C  =  constant  function,  entire 
o 


?1  =  (D  +hPqB)  .(u,PqA  +C) 


3Pi  -i  3P„ 

sr-<»+“p„B>  -<poA+“  ifu 


1  or  , 

+  (D+HPB)'  (P  B+ti  B)(D+tlP  B)"  *(nPA+C) 
O  O  OW  o  o 


The  analyticity  of  PQVP.  implies  analyticity  of  P^  if 


0  £  n  < 


||d"1pob!! 


(2.26) 


Now  assume  exists  and  is  continuous  for 


Then , 


0  ^  u  ^ 


i|D  Pk_1Bl 


—  =(D+,PkB)  -<PkA+*  ^A) 


-1  -1 

+  (D+uPkB)  .(PkB+H^B)(D+uPkB)  (uPkA +C) 


which  exists  and  is  continuous  if 


I 


64 


OS|i< 


From  Chapter  1,  if  condition  (1.21)  is  satisfied,  then  P^ 
converges  to  the  dichotomic  solution.  Thus,  there  exists  a  constant  C 
such  that 

!jp  ||  Vk  s  o 

Thus ,  for 


0  S  n< 


(2.27) 


is  a  holomorphic  function  of  y.  Vk  S  0 

Thus,  P  possesses  a  unique  convergent  MacLauren  series  given 
tc 

by  (2.25)  for  every  u.  bounded  by  (2.27)  Vk  *  0. 

In  a  completely  analogus  manner,  the  series  solution  to 

Vl  *  \  *  (2-28) 

p  =  ubd”1 

o 


can  be  obtained  in  the  form 

P  »  P°  +  HP1  + .  (2.29) 

where 

P°  =  0  (2.30) 

P1  =  BD'1  (2.31) 


pN 


(-  Z  PN'jCPj)-  d'1+apn‘1d'1 

j-1 


(2.32) 


o> 


and  likewise,  (2.29)  is  convergent  for 


0  £  m,  <  - - 

IId'1!!  !lcl! 


a 


Vk  s  0, 


since 


0(n). 


o> 
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C.  Fundamental  Sets  of  Solutions 

In  this  section  we  develop  some  basic  properties  relating  the 
dichotomic  dominant  left  and  right  eigenspace  transformations  of  Chapter  1. 
The  need  for  these  properties  will  become  apparent  in  the  next  section. 

One  of  the  basic  properties  of  linear  homogeneous  systems  of 
differential  equations  of  the  form 

x  *  Fx  (2 .34) 

is  that  a  fundamental  set  of  solutions  [21]  is  of  the  form 


X(t) 


Ft 

=  e 


(2.35) 


For  a  given  initial  value  problem  x(tQ)  =  xqJ  the  solution 
to  (2.34)  for  t  2  tQ  is  uniquely  given  by 

x(t)  *  X(t)X(t  )'1x  (2.36) 

o  o 

or  when  using  (2.35) 


F(t-to) 

x(t)  =  e  x 

o 


Another  property  of  the  fundamental  matrix  (2.35)  is  that 
given  any  nonsingular  matrix  M, 

Y(t)  =  eFt*  M  (2.37) 

is  also  a  fundamental  matrix. 

For  system  (1.1)  satisfying  condition  (1.21)  we  have 

A  ^ 

established  transformation  matrices  P,  Q,  P,  and  Q  such  that 
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Thus , 


— 

A 

B 

= 

C 

D 

— 

— 

_ 

— 

A 

B 

3 

C 

D 

I  Q 
-P  I-PQ 

I-PQ  P 

-Q  I 


A-BP  0 

0  D+PB 

A-PC  0 

0  D+CP 


I-QP  -Q 
P  I 

I  -P 

A  A  A 

Q  I-QP 


A  B 
C  D 


I  Q 

-P  I-PQ 

is  a  fundamental  matrix  for  (1.1). 
Thus,  by  (2.37),  so  is 


(A-BD)t 


(D+PB)  t 


I-QP  -Q 


(2.38) 


(2.39) 


(2.40) 


A  B 
C  D 


or,  the  columns  of 


r 

O' 

Ll_ 

f 

O' 

lT 

e(A-BP)t  Q  ~ 

= 

-P  I -PQ 

_  _j 

.  (D+PB ) t 

-P  I-PQ 

0  e 

(2.41) 


X(t)  = 


e (A-BP)t  Qe(D+PB)t 


.pe(A-BP)t  (I-pg)e(D+PB)t 


(2.42) 


form  a  fundamental  set  for  the  system  (1.1).  Likewise,  using  a  similar 
argument  for  (2.39),  the  columns  of 


X(t)  = 


(I-PQ)e(A-rc)t  Pe^P)t  I 


(2.43) 


-Qe(A-PC)t 


(D+CP)t 


Now,  by  the  dichotomic  property  of  the  transformations,  there 


exist  nonsingular  matrices  T^,  T^,  and 

T.  such  that 

4 

T^1(A-BP)T1  = 

A1 

(2.44) 

t“1(d+pb)t2  = 

A2 

(2.45) 

T^A-PC)!^  = 

A1 

(2.46) 

t^1(d+cp)t4  = 

A2 

(2.47) 

where  is  the  dominant  eigenvalue  matrix  and  A^  is  the  eigenvalue 
matrix  consisting  of  the  rest  of  the  spectrum  of  (1.1). 

Using  (2.44),  (2.45),  (2.46),  and  (2.47)  in  (2.38)  and  (2.39) 


gives 


A  B~j 

A  A 

(I  PQ)TX 

A 

PT 

ri2 

>1 

►— * 

o 

_ 1 

T~^  -T~^P 

L1  1 

!  = 

c  dJ 

_  ^T1 

T2. 

2 

T2^Q  t'^I-QP) 

(2.48) 


A  B 

t3  QT4 

r 

o 

H 

T^tt-QP) 

-t3xq 

-pt3  (i-pq)t4 

-1 

-1 

C  D 

0  A2 

T,  P 

_  4 

T4  __ 

(2.49) 


which,  by  definition  identifies 


(I-PQ)T1 

A 

^2 

T3 

1 

H* 

O' 

-QT1 

T2_ 

* 

_-pt3 

(I-PQ)T4 

(2.50) 


as  eigenvector  matrices  for  (1.1). 
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While  the  magnitudes  of  eigenvectors  are  not  unique,  there 


directions  are.  Thus, 


^2  A2t  T3  Ajt 

e  1  ,  e  1 


(2.51) 


also  serve  as  a  fundamental  set  for  (1.1),  or  in  matrix  form 


y  „  y 

T3e  PT2e 

y  y 

Fr3e  T2e 


(2.52) 


However,  postmultiplying  (2.52)  by  the  nonsingular  matrix 


T3  0 


0  T„ 


(2.53) 


gives  (2.52)  as 


Y(t)  = 


,(A-BP)t  pe  (D+CP)  t 


Pe  (A  BP)  (D+CP)t 


Also,  by  a  similar  argument  on  (2.48)  and  (2.49), 


(2.54) 


(I-PQ)e 


(A-PC)t 


(D+PB )  t 


-<5e(A-k)t 


(2.55) 


also  qualifies  as  a  fundamantal  matrix  for  (1.1). 


Thus,  in  this  section  we  have  established  the  existence  of 
four  fundamental  matrices  for  (1.1)  based  on  the  dichotomic  transformation 
matrices  P,  Q,  P  and  This  flexibility  will  prove  valuable  in  the  next 
section  concerning  asymptotic  expansions  of  our  singularly  perturbed 
model  (2.3). 
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D.  Solution  by  Asymptotic  Expansion  Via  the  Method  of  Vasil 'eva 

Using  the  results  established  in  the  first  two  sections  of  this 
chapter,  we  will  attempt  to  solve  (2.3)  using  asymptotic  expansion  techniques. 

In  [20] ,  the  method  of  matched  asymptotic  expansions  was  proposed 
as  a  method  of  obtaining  asumptotic  solutions  to  the  general  nonlinear 
singularly  perturbed  initial  value  problem 


=  f(z,y,t)  y (0)  =  yQ 
^  dt  =  F^z,y,t)  z(°)  =  z0 


(2.56) 


To  use  this  method,  it  is  assumed  that  the  root  z  =  cp(y,t) 
of  the  equation 


F(y,z,t)  =  0 


(2.57) 


is  stable  in  the  first  approximation  or  specifically,  the  real  parts 
of  the  roots  of  the  characteristic  equation 


DET(  |f 


\I)  =  0 


z^Cy.t) 


(2.58) 


be  negative  in  D,  where  D  is  a  closed  bounded  domain  in  the  variables 
tQ  s  t  ^  t^,  j|z||  <  K^,  |iy|!  <  and  0  £  u  <  Under  this  assumption, 
the  method  can  be  applied  to  (2.3)  as  clearly  carried  out  in  [20]. 

In  our  case  (2.56)  reduces  to 


kL 

dt 


=  Ag  +  Bz 


dz  _ 
dt  " 


Cy  +  Dz 


y(0)  =  yQ 
z(0)  =  z 

o 


(2.59) 


and  Che  assumption  becomes 


Re(X^)  <  0  VX^so(D)  i=l,  ...  m 

with  this  assumption  satisfied,  we  can  proceed  with  the  asumptotic  solution 
method  proposed  by  Vasil' eva. 

A  solution  to  (2,59)  is  sought  in  the  form 


y  =  y  +  tt 

y 

Z  =  Z  +  tt 

z 

where 

y  =  yQ(t)  +  M-y1(t)  +  . 


(2.60) 

(2.61) 


(2.62) 


denotes  a  formal  power  series  in  p.  whose  coefficients  depend  on  t, 
and 

TTy  =  TrQy(T)  +  nn^y  (t  )  + .  (2.63) 


Denotes  a  formal  power  series  in  u>  whose  coefficients  depend  on  p.T  =  t / u, - 
Substitution  of  (2.62),  (2.63)  and  the  analog  bus  expansions 
into  (2.59)  yields 


^  df  +  dT  ^  =  uA(y+rry)  +  w.B(z+tt  ) 
u.  +  TT  ttz  =  C(y+rry)  +  d(z+ttz) 


(2.64) 


Equating  the  coefficients  of  equal  powers  of  p.,  those  depending  on  t  and 
those  depending  on  t  being  treated  separately,  we  arrive  at  the  following 
equations  for  the  variations . 


I 

■A* 


Zeroth  order, 


First  order, 


order 


CyQ(t)  +  DzQ(t)  =  0 


•rjff  z  =  Crry  +  DiTz 
d"  o  o  o 


dy0(t) 

“dt - Ayo(t)  +  B2o(t) 


—  TToy(T)  =  0 


dz.(t) 

+  “]<■> 


drr1  z  (t  ) 

— -  -  CTTjyOr)  +  Dtt^z(t ) 


dy/t) 

-  =  Ay1(t)  +  Bz1(t) 


dn  y(T) 

— -  *  ATToy(T)  +  BttoZ(t) 


dzk-l(t) 

— -  Cyk(t)  +  D2k(t) 


drr,  z(T) 

17 - -  CV(T)  +  DV(T) 


dyk(t) 

“dt  =  Ayk(t)  +  B2k(t) 


(2.65) 


(2.66) 


(2.67) 


(2.68) 


(2.69) 


(2.70) 


(2.71) 


(2.72) 


(2.73) 


(2.74) 


(2.75) 
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dtTky(T) 

dr 


Arrk_ly(T)  +  BTTk-iz(T) 


(2.76) 


Since  we  are  considering  the  initial  value  problem,  the 
matching  conditions  become 


tt  z(0)  +  z  (0)  =  z° 
o  o 

tt  y(0)  +  y  (0)  =  y° 
0  o 

and  for  k  s  l 

rrkz(0)  +  zfc(0)  =  0 
rrky(0)  +  yk(0)  =  0 


(2.77) 

(2.78) 


(2.79) 

(2.80) 


and,  due  to  our  stability  assumption, 

nky(®)  =  rrkz(®)  =  0  k  a  0  (2.81) 


Solutions  of  this  type  are  referred  to  as  "inner"  and  "outer", 
"fast"  and  "slow",  or  "steady  state"  and  "boundary  layer"  depending  on 
the  author. 

Our  purpose  here  is  to  show  that  the  series  solution 


y  =  y  +  rry 

Z  =  Z  +  TTZ 


(2.82) 


are  equivalent  to  the  solution  of  (2.3)  obtained  using  the  fundamental 
matrix  (2.54).  In  other  words, 


] 

?! 
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e(A-BP)t 


(2.83) 


(D-»CP)t 
e  c„ 


(2.84) 


First,  we  seek  an  asymptotic  solution  to  yg^ow  of  the  form 


x(t)  =  xQ(t)  +  M*x(t)  + 


(2.85) 


and  a  solution  of  z  ,  in  the  form 
slow 


-Px(t)  =  -P0xo(t)  -  u(P1xo+P0x1) 


where  x(t)  is  the  transformation  variable  of  (1.82).  When  (2.3)  is 
used  as  the  system  model 


A-P(mOB 


-  +  BP  00 

U 


(2.86) 


Substituting  in  the  formal  power  series  (2.85)  and  (2.7) 
into  (2.86)  and  equating  like  powers  of  u,  we  obtain  the  following 
equations  of  the  variations 


x  *  (A-BP  )x 
o  o  o 


(2.87) 


*1  =  (A-Bpo)xrBPiXo 


(2.88) 


xk  -  <A-BWB  pk-jxj 


(2.89) 
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We  will  now  show  that  the  differential  equations  necessary  to 
solve  for  ,  k  2  0  are  equivalent  to  (2.89)  Vk.  The  equivalence  of 
z  and  (-Px)  is  a  byproduct  of  the  derivation  for  the  equivalence  of  y 

K.  K  K 

and  x,  . 
k 

From  (2.65)  and  (2.67) 


z  =  -  D  *Cy 
o  J  o 


(2.90) 


=  -  P  y 
oJ  o 


dT  ■  (A-BPo)yo 


(2.91) 


Now,  from  (2.69)  and  (2.71) 


_  .  dz 

71 =  Vi  +  D'  df 


(2.92) 


,  dy 

-•  Vi-^df 


—  -2  -l  — 

=  -  p  y  -  D  C (A - BD  C)y 
ol  o 


‘  Vl  -  Vo 


3 r  ■  <A-6Po)yi  -  BPiyo 


(2.93) 


and  thus,  for  (2.73)  and  (2.75) 


z  M  -  p  y  +d 

k  ok 


•i  dVi 


dVi  -  .  dyk-i  r  dV 

dt  o  dt  1  dt 


k-1  dt 


-7T  (2-94) 


z  =  -  P  y 
k  cr  k 


-  d^p  a  +  d_1p  bp 

O  0  0 


- l  -i  k"* 

-D  P.  .A+D  Z  P.BP,  ,  . 

j=Q  J  k-l-j 


•  Vk  •  Plyk-1 

k 

-  ^  P  y. 

j-o  k-J  J 


-  P.y 
k  o 


dt 


(A-BPo)yk  -  B  ^  Pk_.y. 


k-1 

Z 

j=0 


77 


(2.95) 

(2.96) 


which  is  equivalent  to  (2.89)  Vk.  Plus,  it  is  obvious  from  (2.95)  that 

\  =  -  (P*)k> 

We  now  will  show  (2. 84). 

Using  the  dominant  right  eigenspace  iterations,  the  singularly 
perturbed  model  (2.3)  is  transformed  into 
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P  is  obtained  from  (2.78)  and  Q  is  then  obtained  from  (2.90). 
The  dichotomic  nature  of  our  fast  and  slow  components  led  us  to  identify 
fast  and  slow  components  of  y  and  z  as 


^s low  *  <r**>* 

y ,  .  *  Pai 

fast 

zslow  * 

z _  3  'X 

fast 

The  differential  equation 

for  the  fast  state  vector  is 

n+cp  .  . . 

tu  = 

u 

(2.98) 

Let  t  =  t/u,  then  (2.98)  becomes 


^  =  (D+CP)uj(t) 
d^ 


(2.99) 


We  now  seek  an  asymptotic  solution  to  zfast  of  the  form 


'u(T)  =  03  (T)  +  Ud)  (T)  + 
o  1 


(2.100) 


and  a  solution  to  y,  .  of  the  form 
J  fast 


Ru(T)  =  P^Cr)  +  ^(P^Ct)  +  P1«jo(t))  + 


Substituting  the  formal  power  series  (2.100)  and  into  (2.29) 
and  equating  like  powers  of  u,  we  obtain  the  following  equations  of 
the  variations 


d£  (t) 

— - -  =  Dui  (t) 

dT  ov 


(2.101) 


We  will  now  show  that  the  differential  equations  necessary 


to  solve  for  tt^z(t),  k  £  0  are  equivalent  to  (2.103)  k.  The  equivalence 
of  tt,  y (t)  and  ( Ptu )  is  a  byproduct  of  this  derivation. 

K  iC 

From  (2.68)  and  (2.81) 


Thus , 


tt  y  =  0 
o 


dTT  Z 
_ 0 

dT 


—  =  Drr  z 


(2.104) 

(2.105) 


From  (2.72) 

drr  y(T) 

— -  =  Att  y  +  Btt  z  (2.106) 

dT  o  o 


Thus, 


TTxy(T)  =  77^(0)  +  J  [ATTQy(a)  +  Bt7^z(ct)  Jdcr  (2.107) 


To  establish  rr^y(0). 


0  =  77^(0)  +  J  [ATTQy(a)  +  Btt qz(ct) 


(2.108) 


Thus , 


TT^C")  *  -  J  [ATToy(cr)  +  BTToz(CT)]da 


(2.109) 
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Since 


and  as  a 


Finally, 


using 


nQy(a) 


=  0 


TTiy(T)  =  '  I  BTToZ(CT)dCT 


(2.110) 


=  -  f  bd 


,  dn  z(a)da 

■1  O _ 

da 


result 


BD-1[tT  z(“)  -  TT  Z  (t)] 
o  o 


-1 

BD  TT  z(t)  =  P.TT  z(t) 

O  10 


(2.111) 


dn  z(t) 

— — =  DV(T)  +  cpiV(T) 


th 


for  the  k  variation 

CO 

TTky(T)  =  -  P  [Att^_  Ly  (a )  +  Bnk_1a(CT)]dCT 


(2.112) 


(2.113) 


,  00  dn,  ,z(a)da 

■  -  »l  J 


-  J  (A-BD“1D)P1TTk_2z(a)da 


-  J  (A-BD"1C)P2nk_3z(a)da  -  ....  -  J  (A-BD^OP^^z^da  (2.114) 


-1  dTTk-lz(T)  -i  k“2  . 

Vi<T)=D  "dP  c^Vi-jY^) 


(2.115) 


81 


Tr,y(T)  = 

t C 


-  p. 


f*  dnk-iz(o> 


da 


da  -  P„ 


dTTk.22(<J) 


da 


da  - 


“  drr  z(a) 

A  f  O 

-  \  l 


=  +  P2rrk-2z(T)  +  + 


1  k-1 
k-i 

E  P  .tt.z(t) 


k  o 


j=0 


k-j  j 


Thus , 


dTTkz(T) 


k-1 


dr  tmk.(T,  +c  Pk..n.z(T) 


Which  is  equivalent  to  (2.103)  Vk  ^  0.  Also,  from  (2.117), 


(ny)k  -  (Pu)k  Vk  £  0 


Thus,  we  have  shown  that  y  satisfies 

o 


and  that 


dy  _ 

dT  ’  <A  -  BP)po 


dn  (T) 

—j - (D  +CP)nz(T) 


likewise 


TT  =  PTT  (T) 

y  2 


Zn  =  ' 

0  O 


(2.116) 


(2.117) 


(2.118) 


(2.119) 


(2.120) 


(2.121) 


(2.122) 
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The  matching  conditions  (2.77-80)  thus  reduce  to 


yQ  *  y(0)  +  ny(0)  (2.123) 

z  =>  z(0)  +  ttz (0)  (2.124) 

o 


However,  from  (2.117)  and  2.95) 

y0  =  y(°)  +  Ptt2(0) 

zQ  ~  -  Py(0)  +  ffz(O) 


M 

A 

I  p 

y(0) 

Z 

= 

-p  I 

TT  (0) 

0 

_ 

Z 

Thus, 


_  _ 

-1 

— * 

y(0) 

= 

i  p 

y 

■'o 

_V0)_ 

-p  i 

l 

z 

L-  o  — 

(2.125) 

(2.126) 


(2.127) 


(2.128) 


Which  uniquely  determines  the  solutions  to  (1.20)  and  (1.21)  and  then 
(2.83)  and  (2.84).  However,  given  the  fundamental  set  (2.54)  and  the 
initial  conditions  y(0)  =  yQ,  z(0)  =  zq,  the  solution  to  the  initial 
value  problem  is  uniquely  determined  by 


_  _ 

y(t) 

e(A-BP)t  pe (D+PB)t 

-1 

^0 

z(t) 

_pe(A-BP)t  e(D4€P)t 

"P  I 

z 

—  J 

_  _ 

-  — 

_  o _ 

(2.129) 


In  conclusion,  given  the  initial  value  problem  (2.3)  that 
satisfies  condition  (1.21),  the  presence  of  the  singular  perturbation 
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parameter  u  suggests  seeking  a  series  solution  in  two-time  scales.  The 
time  scale  t  corresponds  to  the  slow  or  "outer"  solutions.  The  time  scale 
T  =  t/M-  corresponds  to  the  fast  or  "inner"  solution  which  effectively 
blows  up  the  initial  region  of  nonuniforaity .  Vasil'eva's  method  of 
matched  asymptotic  expansions  is  an  asymptotic  method  of  decomposing  the 
solution  vector  of  (2.3)  into  fast  and  slow  components.  The  mechanism 
for  separation  is  based  on  the  functions  f  and  F  dependence  on  the 
independent  variables  t  and  t .  What  we  have  shown  here  is  that  for  linear 
time -invariant  singularly  perturbed  systems,  the  explicit  time-scale 
decomposition  of  matched  asymptotic  expansions  is  equivalent  to  the 
e i gens t rue ture  decompositions  of  Chapter  1. 
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E.  Simplified  Iterative  Schemes 

One  of  Che  computational  drawbacks  of  the  dominant  left  and 
right  eigenspace  iterations  is  the  computation  of  the  inverses 

(D+uP^)’1  (2.130) 

(D - CP^)*1  (2.131) 

at  every  iteration. 

Looking  at  the  case  (2.6),  the  iterative  matrix  recursion  is 

Pk+1  =  (D+uPkB)'a  •  (C+U-PkA)  (2.132) 

which  can  be  expressed  as 

Dpk+i "  •  “V  Vi  +  c  + 11 V  <2-133) 

if  this  is  approximated  by 

DPk+1  =  -  uPkBPk  +  C  +  U.PkA  (2.134) 

Pk+1  -  D"  l(  C  +4?kA  -  U?kBPk)  (2 . 135) 

will  have  eliminated  the  need  for  the  inverse  (2.130)  at  every 
iteration.  However,  the  question  remains  as  to  whether  (1.135)  converges 
or  not  and  if  it  does,  does  it  converge  to  the  dichotomic  solution  P. 

To  answer  these  questions,  let  us  take  an  in  depth  look  at  the  "equations 
of  variation"  for  (1.132).  In  Section  B,  we  obtained  an  asymptotic 
expansion  in  M-  to  the  .solution  of  (1.132).  We  did  this  by  equating  like 
powers  of  M-  in  the  equation 
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[D+KP°+up£  + .  )B)(P°+1 +BpJ+1  + 


+  +(p£+*p£  + 


) 

)A 


(2.136; 


and  obtaining  the  equations  of  variation 


DPk+l  "  C 


P°  =  D_1C 


(2.137) 


OP*  +  P?BP°  =  P?A 
k+1  k  k+1  k 


opk+i +  ?l  Epi« +  pX«  ■  =kA 


p2  =  o 
o 


Dpf  +  E  P^3pjJ"1'j  =  P^_1A 

k+1  ._n  k  k+1  k 

j  =0 


(2.133) 


if  we  plot  the  solutions  of  the  equations  (2.133),  we  obtain  the  results 
of  Figure  2.1.  In  other  words,  the  variation  does  not  reach  its 
equilibrium  value  until  k  =  N. 
thus,  KN  s  0 


where 


N  .  .  „ 

P  *  £  PJM-J  +  O(^) 

j=0 


(2.139) 


P°  -  d'*c 


_1  -l_oA  -1_0 ,  „0 

P  =  D  PA-D  P  BP 


pN  ,  D-lpN-lA  .  D-1  N‘*  pjBpN-l-j 

j*o 


(2 . 140) 
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In  Section  B  of  this  chapter  we  established  bounds  on  u  such  that 
(2.139)  is  a  convergent  MacLauren  series. 

Now,  let  us  obtain  an  asymptotic  scales  solution  to  (2.135)  and 
compare  it  to  (2.139).  Substituting  in  a  formal  series  for  P^,  we  obtain 

<?k+l  +^+l  +  ••••>"  D^tc  •  •  -)A 

-  ii(P°+U?J  +  ...)B(P°+^pJ +  ...)]  (2.141) 


which  gives  us  the  following  variational  equations 

p£+1  =  D_1C  (2.142) 

— 1  “1—0  “1-0  ~o 

P.  V,  *  D  lPTA  -  D  X  BP, 
k+1  k  k  k 


.  =  d“1b£“1A  -D_1  S  P^BP^-1^ 
k+1  k  j=0  k  k 


(2.143) 


If  we  plot  the  solutions  of  these  equations  we  obtain  the  results  of 


Figure  2.2. 
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Again,  the  N  variation  does  not  reach  its  equilibrium  value  until  K  =  N. 
Moveover,  Vu  2  0 

~  ^  4  j  2 

P  -  S  p-V  +  0(u  )  (2.144) 

j=0 

where 


p°  =  d'1C 


~1  -1  O  1~0  ~0 

P  =  D  P  A  -D  p  BP 


F  =  -  D'1  2  PjBPN_1"j  (2.145) 

j=0 

which  are  equivalent  to  (2.139)  2  0.  However,  we  know  that  if 

0<K  - i-i -  (2.146) 

1Id~  1 II  ItBii  M 

Ihen  (2.139)  is  a  convergent  MacLauren  series.  Thus,  for  u  satisfying 
(2.146),  (2.143)  is  bounded  and  converges  to  the  dichotomic  solution  P. 

Obtaining  bounds  for  the  convergence  of  (2.135)  to  the 
dichtomic  equilibrium  were  first  established  in  [4].  In  this  case  a 
contraction  mapping  argument  is  used  to  show  convergence  if 


where  A  =  A 
o 


BD_1C 


_ 1 _ 

3i!D“1||(i!Aoll  +j|Bi|  lb_1cil 


0  £  h 


(2.147) 
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To  show  one  improvement  over  this  bound,  we  note  that  from  [4] 


,  2  iiA  ii  |:D"1cl! 

•  Pi!  s  iiD_1c||  + - 2 - 


Ia  i!  +  ;Ib||  llD_1ci! 

o 


(2.148) 


l„n2 


3iIa  il  lio"  1c|!  +  i'bIi  I!d_1c! 

s  - 2 - 

IIaJI  +l|Bii  i,D"1c|i 


(2.149) 


thus  bound  (2.146)  becomes 

1 


0  ^  u 


IIaJI  +  i f B i j  !|D_1ci! 


-1 


If  ijs !■  (3 [Ia  II  I|d'1c||  +Hb||  IId"1^2) 


||a  II  +  j | b I !  !Id‘1c| 

o _ ' 


d"1!!  IIbII  !lD'ic!|(3|k  II  +||b||  d"xc||) 


-i 


-l 


(2.150) 


which  is  obviously  less  conservative  than  (2.147).  The  dual  to  this 
result  involving  equation 

Vl  =  li(B+AV(D+C\)"1  (2.151) 

is  derived  in  the  analogus  manner.  The  simplified  interation  in  this 
case  becomes 

Pk+1  =  (UB  +  uAPk  -  PkCPk)D"1  (2 . 152) 

P  =  U.BD-1 
o 

and  converge  to  the  dichotomic  solution  P  when  p  is  bounded  by  (2.33). 

In  the  next  section  we  give  an  example  highlighting  many  of  the 
important  results  in  this  chapter. 
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F.  Example  -  Fundamental  Solution  Sets 

Given  the  same  power  system  model  from  our  previous  example 
given  here  in  singularly  perturbed  form 


r  ~ 

— 

— 

— 

-.2 

0 

0 

0 

0 

0 

0 

0 

| 

0 

-.167 

0 

0 

0 

0 

0 

.167 

:: 

! 

y 

0 

0 

-.5 

0 

.2 

0 

0 

0 

.  y 

t 

0 

.009 

0 

-.112 

-.063 

.014 

.116 

.01 

0 

-.075 

10.0 

-9.101 

-3.994 

-.112 

-.927 

-.08 

0 

2.0 

0 

0 

0 

-2.0 

0 

0 

; 

az 

0 

0 

0 

-.278 

1.319 

0 

-1.389 

0 

z 

i 

4.75 

0 

0 

0 

0 

0 

0 

-5.0 

i 

(2.153) 


We  seek  a  solution  of  the  form 

(A-BP)t  -  (D+CP)t 


Pe 


-Pe 


(A-BP)t 


(D  +CP)t 


► 

— 

— 

-1 

— 

i 

* 

p 

• 

. 

X 

0  ! 

-p 

i 

i 

i 

t 

_ 

! 

(2.154) 


using  the  reduced  order  iterations  of  Section  E. 

To  use  (2.135)  and  (2.154)  we  must  first  examine  the  bounds 
on  u  that  guarantee  convergence . 

For  (2.135), 


0  £  u  <  3.56285 


and  for  (2.152) 


0  £  u  <  4.42141 

Since  u  =  1  in  our  problem  (2.153),  we  are  well  within 


(2.155) 


(2.156) 


bounds . 
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Using  ? 


A  ■  P  B  = 
o 


-0.20000 

0.15834 

-0.00312 

0.00877 


0.00000 

-0.16667 

-0.00766 

0.02153 


0.00000 

0.00000 

-0.08981 

0.09635 


0.00000n 

0.00000 

-0.36571 

-0.22145 


(2.157) 


with  eigenvalues 


-0.15563  +0.17579  J 
-0.15563  -9.17579  J 
-0.16667  +0.00000  J 
-0.20000  +0.00000  J 


(2.158) 


Likewise,  using  P 


D+PC3 

o 


-4.31691 

0.00000 

1.32208 

0.00000 


-0.03023 

-2.00000 

0.00179 

0.00000 


0.04757 

0.00000 

-1.36749 

0.00000 


-0.05415 

-0.06667 

0.00051 

-5.00000 


(2.159) 


with  eigenvalues 


-4.33808  +0.00000  j 
-1.34632  +0.00000  J 
-2.00000  +0.00000  J 
-5.00000  +0.00000  j 


(2.160) 


Using  an  initial  condition  of 


xC  33  (1,  2,  3,  -  2,  1,  -1,  4,  2) 
o 

The  full  accuracy  original  states  of  (2.153)  are  plotted 
versus  their  Oth  order  approximations  using  (2.145)  on  the  next  several 
pages.  Some  higher  order  plots  are  also  given  where  convergence 
verification  is  needed. 


Actual  State 
Approximated 


</}f-  Cl-lii 


Figure  2.5.  State  3  0 


csoNooas: 


10. 


Figure  2.10.  State  5  1 


Figure  2.13.  State  7  1  order  approximation. 
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D  +  CP  = 
o 


-1.65397 

0.00000 

0.44248 

0.00000 


0.04443 

-0.66667 

0.00179 

0.00000 


0.66609 

0.00000 

-0.44156 

0.00000 


-0.00082 

-0.06667 

0.00051 

-1.66667 


(2.166) 


with  eigenvalues 


-1.86153  +0.00000  J 
-0.23400  +0.00000  J 
-0.66667  +0.00000  j 
-1.66667  +0.00000  J 


(2.167) 


Some  eigenvalues  are  very  accurate  while  some  have  a  sub¬ 


stantial  error.  Using 


3 

“0.20000 

0.00000 

0.00000 

0.00000 

0.17993 

-0.16667 

0.00000 

0.00000 

A  -  BP5  * 

-0.00090 

-0.00270 

-0.05903 

-0.32775 

with  eigenvalues 

-0.00053 

0.02545 

0.11841 

-0.02838 

-0.04371  +0.19641  J 
-0.04371  -0.19641  J 

-0.16667  +0.00000  J 
-0.20000  +0.00000  J 


(2.169) 


likewise,  using  P^ 


D  +  CP5  = 


-1.90740 

0.00000 

0.44430 

0.00000 


0.00077 

-0.66667 

0.00253 

0.00000 


-0.11943 

0.00000 

-0.43118 

0.00000 


-0.00521 

-0.07408 

0.00040 

-1.66667 


(2.170) 


with  eigenvalues 


-1.87053  +0.00000  J 
-0.46805  +0.00000  J 
-0.66667  +0.00000  J 
-1.66667  +0.00000  J 


(2.171) 


Continuing,  the  eigenvalues  of  A  -  BP^  are 
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-0.03291  +0.21694  J 

-0.03291  -0.21694  J 

-0.16667  +0.00000  J 

-0.20000  +0.00000  J 


and  the  eigenvalues  of  D  +  CP^  are 


-1.88648  +0.00000  J 
-0.45457  +0.00000  J 
-0.66667  +0.00000  J 
-1.66667  +0.00000  J 


Thus,  we  have  convergence.  However,  the  rate  of  convergence  decreases  as 
approaches  the  circle  of  convergence. 


.172) 
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III.  CONCLUSIONS 

In  this  report  we  have  presented  various  methods  used  decomposing 
large  scale  systems  into  reduced  order  subsystems.  These  methods  came  under 
such  names  as  quasi  steady  state,  Riccati-iterates ,  and  matched  asymptotic 
expansions  in  previous  works.  Using  power  iterations  for  obtaining  the  left 
and  right  dominant  eigenspace  of  a  matrix  we  have  been  able  to  show  that  the 
convergent  behavior  of  the  above  methods  to  be  equivalent.  Decomposition  of 
time  scales  in  linear  systems  reduces  to  separation  of  invariant  eigenspaces. 

If  there  is  a  magnitude  separation  in  eigenvalues,  then  it  is  possible,  using 
the  above  mentioned  procedures ,  to  iteratively  solve  for  the  invariant  eigen¬ 
spaces  and  generate  reduced  order  models. 

In  Chapter  II  of  this  report  we  have  shown  that  a  linear  singularly 
perturbed  system  that  is  decomposed  through  matched  asymptotic  expansions  is 
in  essence  an  eigenspace  decomposition.  The  presence  of  the  singular  pertur¬ 
bation  parameter  u  enables  us  to  obtain  convergent  power  series  solutions  to 
the  block  diagonalization  matrices  P  and  P  that  enabled  us  to  obtain  computa¬ 
tionally  more  efficient  algorithms  for  obtaining  the  time  scale  decompositions. 
The  singularly  perturbed  model  can  thus  be  looked  upon  as  an  application  of 
parameter  imbedding  in  an  effort  to  obtain  a  power  series  decomposition  in  two 
time  scales,  namely  t  and  t«  t/p. 

Finally,  in  the  Appendix,  we  give  an  application  of  this  work 
involving  partial  pole  placement. 


APPENDIX  -  APPLICATIONS  TO  PARTIAL  POLE  PLACEMENT 


There  are  many  applications  using  the  techniques  developed  in  the 
first  two  chapters.  They  include  robust  designs,  reduced  order  regulator 
designs,  and  reduced  order  modeling  only  to  mention  a  few.  In  this  appendix 
we  will  show  the  time  scale  decomposition  techniques  can  be  used  to  implement 
partial  or  full  pole  placement  design  on  both  continuous  and  discrete  systems. 

We  will  now  be  considering  the  completely  state  controllable 


system 


u£  RP. 


(Al) 


If  the  open  loop  eigenvalues  satisfy  (1.21),  then  we  can  apply 
transformation  (1.14)  which  transforms  (Al)  into 


y 

"a-bp 

B 

y 

G 

= 

+ 

n 

0 

D+PB 

A 

H+PG 

(A2) 


where  P  is  obtained  using  either  (1.17)  or  (2.135).  The  transformation 
involved  here  can  be  written  as 


(A3) 


which  possesses  the  explicit  inverse 
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Ill 


Observe  now  that  the  pair  (A-PC,  G-PH)  spans  only  the  slow  controllable 
subspace.  Thus,  design  a  feedback  F  such  that 

(A*  +  G*F)  A*  *  k-fc 
G*  -  D-t-CP 

has  N  desired  poles. 

The  control  is  of  the  form 

u  «  F? 

=»  F(y-Pz) 

■  •  -rfl  [j] 

and  the  resulting  closed  loop  system  has  N  eigenvalues  of 

a (A*  +  G*F) 


(All) 


(A12) 


and  M  eigenvalues  of 


o(D  +  CP) , 


In  general,  both  slow  and  fast  modes  may  be  designed  for.  In  this  case  a 
general  two-state  design  procedure  may  be  implemented.  Assume  we  have  used 
either  (1.83)  or  (1.94)  on  (Al)  and  thus  have  the  form 


A* 

0 


X 

1 

o 

* 

+ 

w  _ 

w 

u. 


(A13) 


We  arbitrarily  chose  to  design  for  the  slow  subsystem  first.  Thus,  we  chose 


an  F  such  that 
s 


a(A*  +  G*F  ) 
s 


has  N  desired  "slow"  eigenvalues.  Letting 

u  ■  u  +  u. 
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where  now 

u  ■  F  x, 
s  s 

the  partially  closed  loop  system  looks  like 


- 

— 

"1 

• 

X 

m 

A*  +  G*F 

s 

0 

X 

w 

H*F 

L  s 

* 

o 

w 

Now,  let 


v  *  w  +  Sx 


(A14) 


(A15) 


which  transforms  (A14)  into 


X 

A*  +  G*F  0 

“ 

X 

G* 

a 

s 

+ 

V 

H*F  +s(A*+G*F  )-D*s  D* 

L  s  v  s  -1 

_w 

H*+  SG* 

We  pick  S  such  that 


(A16) 


H*F  +  S(A*+GF  ) -D*S  -  0.  (A17) 

s  s 

This  Lyapunov  type  equation  has  a  unique  solution  if 

a(D*)  <"•  a  (A*  +  GF  )  -  $  (A18) 

s 

and  thus  can  be  solved  algebraically  [12].  If 

inf|c(D*)|  >  sup|a(A*  +  G*F  )|.  (A19) 

s 

Then  the  iterative  scheme  (1.92)  may  be  used  to  solve  (A17) .  With  this  S,  the 
pair  (D*.  H*+SG*)  now  spans  only  the  fast  controllable  subspace,  and  we  can 
design  a  feedback  gain  F^  such  that 

cr(D*+  (H*+SG*)Ff) 


has  M  desired  eigenvalues. 
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Thus,  our  composite  control  is 

u  *  F  x  +  F.v 
s  f 

-  G  x  +  F,(w+  Sx) 

s  f 

*  (F  +  F,s)x  +  F-w 
s  r  r 

-  UFS  +  Ff,):rf][;]  020) 

and  using  either  transformation  (1.84)  or  (1.95),  (A20)  can  be  expressed  in 
terms  of  our  original  state  variables.  This  control  place  N  eigenvalues  of 

o(A*  +  G*F  )  (A21) 

s 

and  M  eigenvalues  of 

cr(D*+(H*+SG*)Ff).  (A22) 

This  technique  has  been  applied  to  singularly  perturbed  systems  [22]  where  it 
is  shown  to  be  a  generalization  to  results  obtained  in  [17,18,19].  This 
technique  is  also  applicable  to  discrete  time  models  as  shown  in  [24],  In  this 
case,  the  dominant  eigenvalues  are  part  of  the  "slow"  spectrum. 


r 
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